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Abstract 

This paper presents our interpretation and understanding of the different 
backgrounds in the EDELWEfSS-I data sets. We analyze in detail the several 
populations observed, which include gammas, alphas, neutrons, thermal sensor 
events and surface events, and try to combine all data sets to provide a coherent 
picture of the nature and localization of the background sources. In light of 
this interpretation, we draw conclusions regarding the background suppression 
scheme for the EDELWEISS-II phase. 
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1 Introduction 



EDELWEISS-I is the first phase of an experiment aiming at the direct detection 
of WIMPs (Weakly Interacting Massive Particles) which could constitute the 
dark matter halo of our Galaxy. It uses cryogenic detectors able to measure 
simultaneously the heat and ionization components of the energy deposit in- 
duced by the elastic scattering of a WIMP off a target nucleus (see e.g. [1] for 
a review). 

The final analysis of the EDELWEISS-I data [2] covers a total fiducial ex- 
posure of 62 kg.d. The limits on the neutralino scattering cross-section were 
obtained from the observation of 40 nuclear recoil candidates with recoil en- 
ergies between 15 and 200 keV, of which 3 are between 30 and 100 keV. The 
limits were obtained without the subtraction of any background, although the 
presence of a coincidence between two detectors and the study of charge col- 
lection distributions suggested that at least some of these events are due to a 
neutron background and surface interactions of electrons. In the preparation 
of the second phase of the experiment, it was necessary to perform a thorough 
investigation of all events in EDELWEISS-I before the nuclear recoil selection, 
in order to better assess their origin and devise means to remove as many back- 
ground components as possible for EDELWEISS-II. 

After a short description of the EDELWEISS-I experiment and of the active 
background rejection capabilities of heat/ionization cryogenic bolometers, this 
paper presents the evaluations of the backgrounds related to different origins : 
gamma radiation, alpha particles and other surface events, and neutrons. While 
the energy range of interest for WIMP detection is limited to below 100 keV, 
high energy gamma and alpha lines can prove very useful to understand the 
backgrounds in the range relevant to WIMP interactions. 

2 EDELWEISS-I 

2.1 General setup 

EDELWEISS is located in the Laboratoire Souterrain de Modane (LSM) which 
provides a ~4800 m.w.e. rock shielding against cosmic muons, reducing the 
vertical muon rate to 4.5 /day/m^. The 320 g germanium detectors are oper- 
ated at a very low temperature (typically 17 mK) within a dilution cryostat. 
This cryostat, made mostly of ultrapure copper, is further protected from ex- 
ternal radioactivity by 10 cm Cu and 15 cm Pb shields [3]. Pure nitrogen gas 
is circulated inside this shielding to reduce radon accumulation. An external 
30 cm layer of paraffin protects the experiment from neutrons created in the 
rock. Inside the cryostat, the detectors are shielded from the radioactivity in 
components of the electronics by 7 cm of roman lead 0- The cold electronic 
components (in total, nine FETs and a dozen of resistors and capacitances) are 
kept as far away as possible from the detectors, and the rest of the acquisition 
chain remains outside the cryostat. Fig. 1 presents a schematic cut view of the 
cryostat inside its Cu and Pb shields. 

^The archeological lead comes from an antique roman ship. The wreck was investigated 
during four campaigns (1984-86) supervised by M. L'Hour of the Direction des Recherches 
Areheologiques Sous-Marines (4,. 
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The detectors themselves are encased into 1 mm thick uhrapure Cu casings, 
and held in place with the help of teflon balls and three small Cu springs (see 
the inset in Fig. 1). The thermal sensor consists of a small (1 mm'^) Neutron 
Transmutation Doped (NTD) germanium crystal glued directly onto each crys- 
tal. The thermal coupling to the cryostat is assured by ultrasonic bonding of 
several micrometric gold wires on gold pads. Electric connections are assured by 
the same kind of gold wires linking the electrodes to contact tracks on the copper 
casing. The wires going up to the 4 K pre-amplifier level are low radioactivity 
coaxial cables. Further details can be found in ^5j. 

2.2 Active background rejection capabilities 

It is impossible to completely shield the detectors from external radiation. 
Most of the radiation reaching the detectors is in the form of gamma rays. In 
EDELWEISS-I, the rate of gamma interactions exceeds the one expected from 
WIMPs by at least a factor 10^. In order to address this, the EDELWEISS 
detector technology offers means to actively discriminate between electron re- 
coils caused by photons and electrons, and nuclear recoils caused by neutrons 
or WIMPs. Indeed, for the same interaction energy, nuclear recoils have an 
ionization efficiency on average three to four times less than electron recoils, 
depending on the energy By measuring simultaneously a heat signal and 
an ionization signal and considering the ratio of both parameters, it is possible 
to reject more than 99.9 % of the gamma interactions while keeping a 90 % 
efficiency for nuclear recoils down to an energy of 15 keV f2l. The value of this 
threshold depends on the experimental resolutions, which have been as low as 

1.0 keV on the ionization channel and 0.3 keV on the heat channel. 

The ionization measurement is made possible by two aluminum electrodes 
(thickness 60-100 nm) deposited onto top and bottom surfaces of the germa- 
nium detector. By applying a moderate voltage between the electrodes, charge 
carriers created by an interaction can be collected and give rise to a signal. The 
selected voltage is 4 V, which is high enough to efficiently collect charges but 
not too high in order to preserve the discrimination capabilities. Furthermore, 
one of the collection electrodes is actually separated into a center part and an 
outer guard ring (see Fig. 2). This allows to define a fiducial volume inside 
the detector where a reliable charge collection is expected, as opposed to the 
lateral sides where electric field lines can escape the crystal before they reach 
the electrodes, leading to an incomplete charge collection. The current NTD 
bolometer technology has been "upgraded" with the addition of a '^60 nm layer 
of amorphous Ge (GGA detector type) or Si (GSA detector type) just below 
the electrodes (see Fig. 2), which essentially gets rid of the majority of surface 
interactions [IllH]. 

3 Gamma background 

3.1 Radioactivity measurements and material selection 

More than 99.5 % of the interaction rate in the EDELWEISS-I detectors is due 
to gamma events. The 15 cm thick lead shield stops most of the gammas from 
outside the experimental setup and the observed gamma background comes from 
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the inside. To select low activity materials a separate low-background counting 
facility was built: a n-type coaxial High Purity Ge diode of 210 cm'^ operated at 
77 K with archeological lead shielding. This dedicated HPGe diode was not set 
up at the time of construction of the shielding and of the cryostat. Thus, only a 
partial material selection could be performed before the experiment; the copper 
used for the shielding had to be measured "a posteriori" . The HPGe diode was 
intensively used for internal material selection close to the detectors. Table 1 
displays the measured radioactivities of some components of the EDELWEISS-I 
set-up, most of them situated in the immediate vicinity of the detectors: copper 
detector holders, copper springs, teflon pads and coaxial cables. 

3.2 Data sets and high-energy data reconstruction 

The present analysis is based on two data sets. The first is the run labeled 2003p 
in Ref. 2i], representing 39.4 kg.d of data (total volume of the three detectors) 
taken in very stable conditions, with similar performances of the three detectors 
in terms of energy threshold and recoil energy resolution. For the WIMP search, 
the amplifier gains, ADC bit ranges and channel sensitivities are optimized for 
the low-energy WIMP signal, with the consequence that ionization signals above 
1 MeV saturate the digitizers. Saturation of the heat signal occurs at 350 to 
600 keV, depending on the detector. This is far above the expected range for 
WIMP signals, but it affects the identification of radioactive backgrounds where 
gamma lines up to 2.6 MeV can prove useful. 

For this reason, a dedicated run, labeled 2003h, has been recorded in the 
same experimental conditions as the run 2003p, but with all amplifier gains 
reduced by a factor 10. The total exposure for this run is 9.3 kg.d in the total 
volume of the three detectors. 

Comparing the data from the runs 2003p and 2003h, it was confirmed that 
the saturated signals of the run 2003p could be corrected and their amplitudes 
could be calibrated reliably. This is possible because, as described in Ref. [2], 
extensive samples of the signal time profiles associated to each event are stored 
onto disk. However, the filtering applied by the acquisition system distorts the 
shape of saturated pulses. Therefore new templates were built by filtering ideal 
events truncated at different levels of saturation. The template used for a given 
event is the one giving the best fit. It is obtained by varying the saturation level 
of the template and minimizing the corresponding ^ . 

This allowed us to reconstruct gamma and alpha events with recoil energies 
as high as 7 MeV, and a full width resolution of about 3 % at 2.6 MeV (ionization 
channel). This reconstruction method suffers however from a lack of sensitivity 
in the case of weakly saturated pulses as the procedure is comparatively more 
sensitive to the baseline level of the events which has a direct impact on the 
saturation level at a given energy. Conservatively, as discussed in [9j , the method 
is not used below 1 MeV for the ionization channel, and below 350-600 keV for 
the heat channel, depending on the detector. 

3.3 Experimental gamma background spectrum 

The energy calibration of ionization signals was performed using ^''Co (122 
keV), i37Cs (662 keV) and ^°Co (1173 and 1332 keV) sources. The full width 



3 



resolution was about 2.5% and the response of the ionization channel was linear 
up to 1332 keV. 

Fig. 3-a shows the ionization speetrum corresponding to the total volume 
of the three detectors for the two data sets 2003h (9.3 kg. day) and 2003p (39.4 
kg. day). In these spectra, the 2003p data above 1 MeV have been reconstructed 
as described in the previous section. As we are concerned with the gamma 
spectrum, the alpha events, which appear with low quenching values essentially 
below 2 MeV ionization energy, have been eliminated by a cut on the quenching 
factor Q{Q < 0.5), see below section 4. The most noticeable feature is a Comp- 
ton backscattering bump around 200 keV. Lines originating from U-Th series 
are visible, notably the ^o^Tl line at 2614 keV with a full width resolution of 
about 3%. A peak of at 1461 keV is also present. 

3.4 Background simulations 

Monte Carlo simulations are performed under GEANT3 [T^ using the geometry 
given in Fig. 1 and the measured activities or limits of Table 1. Several sources of 
gamma background are successively studied and comparisons with the observed 
spectrum are given in Fig. 3-b to 3-d. 

3.4.1 Copper detector holders and thermal shields 

The most massive materials in the immediate vicinity of the detectors are detec- 
tor holders, cryostat structure and thermal shields. They are made of about 20 
kg of ultra-pure copper (OFHC copper, purity > 99.99%). The upper Hmit on 
their U/Th content is 0.1 ppb (Table 1). After machining, the copper holders 
were brought underground only a few months before the start of the exper- 
iment. Therefore, the cosmogenic activation of ^°Co (half-life 5.3 y) at the 
surface has to be considered. The experimental limit on ^"Co is compatible 
with an equilibrium concentration of about 1 mBq/kg [TT]. The U-Th and ^°Co 
contributions, assuming activities equal to the limits, are compared in Fig. 3-b 
with the experimental spectrum. The simulated continuum is too low by one 
order of magnitude, with much more pronounced lines than actually observed. 

3.4.2 Radon 

In the lead-copper shield, near the detectors, about 20 liters of air are trapped 
(see Fig. 1). Although continuously flushed with pure nitrogen emanating from 
the liquid nitrogen dewar, this volume is a potential source of radon contami- 
nation. No radon concentration measurement has been made and the decay of 
■^^"^Rn is simulated assuming an activity of 10 Bq/m"^, the mean value measured 
in the LSM cavity. As can be seen on Fig. 3-c, the simulated continuum is 
again much too low, and the predicted lines are not observed. Even with an 
unrealistic high concentration, the decay of radon can't account for the observed 
spectrum. 

3.4.3 Copper shield 

The inner 10 cm thick copper shield (Fig. 1), is made of about 1 ton of copper 
bricks (electrolytic copper, 99.9% purity), purchased and brought underground 
in the early nineties. Measured "a posteriori", this copper shows a small but 
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measurable contamination in U-Th series (Table 1). The cosmogenic ^°Co is, 
as expected, not detected after about two half-lives of decay underground. The 
simulation shows that the U-Th content of this copper shield accounts for most 
of the features of the background energy spectrum (see Fig. 3-d). Gammas orig- 
inating from the bulk of this copper shield go through several centimeters of 
material before they reach the detectors, and, as a consequence, lines aren't very 
pronounced, Compton diffusion being the dominant energy dissipation process 
above ~ 150 keV. The Compton backscattering bump around 200 keV is re- 
produced though the lower amplitude in the simulation indicates that an other 
source of high energy gamma rays might be present. The rate of the ^oS'pj jjj^g 
at 2614 keV (1.1±0.2 counts/kg/d)is reproduced within 30% by the simulation, 
a satisfying agreement given the uncertainty of 60% on the measured thorium 
concentration in copper. 

3.4.4 

A peak of ^"K at 1461 keV (1.6±0.5 counts/kg/d) is present in the energy 
spectrum. Table 1 shows that the highest *'^K activity is found in the wires 
to the detectors (teflon sheath), but the involved mass is very low (few grams) 
and the simulated rate is two orders of magnitude lower than the observed one. 
The same holds for '*°K contaminations of Cu springs (1.3 g for one detector) 
and teflon balls (0.4 g for one detector) in the detector holders. The present 
measurement of the copper of the gamma shield provides only an upper limit 
(Table 1), corresponding to a rate of 3.5 counts/kg/d in the ^''K peak, which is 
consistent with the measured one (Fig. 3-a). 

3.4.5 Summary 

Radioactivity measurements and associated simulations have shown that most 
of EDELWEISS-I gamma background arises from a tiny U/Th contamination 
of the very massive copper shield rather than from radioactive contaminations 
close to the detectors. 

4 Backgrounds from alpha particles and surface 
events 

4.1 Alpha particles 

One of the most noticeable features revealed in the 2003h data is the presence 
of a distinct population localized at a recoil energy En — 5.33 ± 0.03 MeV 
with a quenching factor H Q = 0.30 ± 0.02, consistent over the three detectors 
(see Fig. 4). This population was later confirmed in the 2003p data using the 
high-energy reconstruction procedure detailed in section 3.2 (see Ref. [5]). The 
rates vary from 2.4 ±0.6 to 5.0 ±0.8 counts/kg/d in the fiducial volumes (center 
electrodes) of the three detectors, and from 13 ± 2 to 25 ± 2 counts/kg/d in the 
lateral volumes (guard electrodes). We observe a significantly higher rate in the 

^The quenching factor is defined as the relative ionization efficiency between a nuclear 
recoil and an electron recoil of the same real recoil energy. 
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top detector than in the bottom detector (Table 2), a fact for which we have no 
explanation. 

Given the energy range and peculiar value of the quenching factor, we link 
this population to the interaction of alpha particles in our detectors. Our expla- 
nation is that the detectors themselves and/or their close environment, i.e. their 
copper holders, suffer from a ^-'^^Pb contamination. This isotope, with a half-life 
of 22.3 years, is a daughter of ^^^Rn. It can be implanted on a surface exposed 
to an atmosphere containing radon during the fabrication and handling of the 
detectors. The last disintegration of the chain ^^°Pb^^-^°Bi^-^^°Po^^°^Pb pro- 
duces an alpha with an energy of 5.3 MeV. If the events with Q < 0.5 on Fig. 4 
are interpreted as alphas, their energy spectrum restricts their origin to either 
the detectors themselves or the copper surfaces (< 1/^m) surrounding them. In- 
deed, given the very low penetration length of alpha particles of such energies 
in germanium or copper (~ 10/im, see Table 3), a volume contamination of any 
material other than the germanium itself would undoubtedly lead to the obser- 
vation of an alpha energy continuum down to MeV. While we do observe such 
a tail to low energies, it is clearly not the dominant feature. 

The localization in quenching is also interesting. Although such a phe- 
nomenon may be associated with incomplete charge collection due to surface 
interactions, it is difficult to explain why we observe a constant value of Q ~ 0.3, 
and not a range of values down from Q = I. Previous studies |141ll5j have shown 
that charge collection efficiency for alpha interactions in Ge is similar to that of 
gammas. However, those measurements were obtained at room or liquid nitro- 
gen temperatures, and under a drifting field of several thousands V/cm. Despite 
the lack of more relevant results, it is not unreasonable to assume that, in our 
case, because of the particular nature of an alpha interaction in germanium 
and the high local density of charges created, some systematic recombination 
before collection takes place. This would typically lead to a constant value of 
the quenching ratio, as is observed in our data. 

4.2 Heavy nuclear recoils 

When ^^"Po decays to ^osp]-, 5 3 ]y[gY alpha is produced and the ^"^Pb nucleus 
recoils with a kinetic energy of 103 keV. Depending on whether this takes place 
at the surface of the detector itself or at the surface of the copper in front 
of the germanium, what we expect to see is substantially different. In both 
cases, the penetration length of such a heavy particle at such a low energy is 
so small (a few tens of nm) that all of the ionization signal is lost, charges 
being absorbed either in the aluminum electrode (100 nm) or in the amorphous 
sub-layer (another 60 nm) if the Pb nucleus hits a surface not covered by an 
electrode. If the contamination is localized on the surface of the detector, the 
heat signal should correspond to the full recoil energy of 103 keV. In contrast, if 
the Pb nucleus originates from the copper, then it has to go through a thickness 
of material equivalent to its implantation depth before it reaches the detector, 
resulting in a partial loss of energy. 

In order to isolate this heavy recoil population, we looked at all the events 
in our data compatible with a signal above detection threshold on the heat 
channel, and below threshold on the ionization channel. This led to a clas- 
sification into three categories: a) sensor "NTD" events, induced inside the 
small thermal sensor by its own radioactivity; b) random noise events, linked 
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for example to anomalous microphonic episodes inside the cryostat, and c) "le- 
gitimate" ionization-less events, among which potentially lie the heavy recoils 
we seek. As described in Ref. [2], NTD events present a very particular pulse 
shape, significantly shorter than regular bulk events. We use shape discrimina- 
tion to eliminate populations (a) and (b) and isolate population (c) (see Ref. [9] 
for details). 

Using the data from run 2003p, in each detector we find a population of 
events clearly contained below 100 keV (Fig. 5), with rates varying from 1.5±0.5 
to 5.4 ± 0.8 counts/kg/d in the total volume of the detectors for recoil energies 
greater than 40 keV (see Table 2). Interestingly, as was the case with alpha 
interactions, the top detector registers more events than the bottom one, and 
in the same ratio within the error limits. This tends to confirm that both 
populations are indeed linked. The fact that we do not see a clear peak at 
103 keV but a roughly uniform front below this energy also implies that the 
contamination is localized exclusively on the copper surfaces, and not on the 
detectors themselves. 

4.3 Surface beta interactions 

If the contamination is indeed linked to an exposure to radon, then we expect to 
see the products of all disintegrations in the ^^°Pb chain, in addition to the 5.3 
MeV alphas and heavy recoils. In particular, the ^^''Pb decay to ^^'^Bi leads to 
a complex spectrum of low energy conversion and Auger electrons together with 
two beta spectra with end-points at 63.5 keV and 17 keV. Due to the value of its 
penetration length in germanium (see Table 3), the 46.5 keV gamma ray does 
not contribute significantly to the surface event budget. The decay of ^-'^^Bi to 
•^^°Po emits another beta electron (branching fraction 100%) with an end-point 
of 1.16 MeV. 

According to Monte-Carlo simulations using the CASINO code [T^ (see 
Table 3), an electron of 100 keV impinging at normal incidence the germanium 
surface will lose 90% of its energy in the first 20 /im. This value goes up to 
~700 //m for a 1 MeV electron. As for the nuclear recoils (sect. 4.2), the 100 
nm aluminum electrode and the 60 nm amorphous semi-conductor sub-layer 
constitute a "dead ionization zone" in our detectors. As for the alphas (sect. 
4.1), there exists a zone under the electrodes where electron-hole pairs are not 
properly collected. The depth of this zone can be as large as 10 /im, as studies 
using detectors with a different design have shown [71 [TB] . 

We therefore expect to observe some incomplete charge collection for a sig- 
nificant part of the events generated by the decay electrons. In the 2003p data, 
we do indeed observe a population of "intermediate" events between the elec- 
tron recoil band and the nuclear recoil band (Fig. 6). As shown in Ref. [2], 
this population is absent when the detectors are exposed to a ^'^^Cs source of 
penetrating gamma rays. In order to quantify this population, we compare our 
low-background data to ^'^^Cs gamma calibration data. This allows us to define 
an area in the (Q, Ei^) plane (with En <80 keV) where we can be confident 
that events are not caused by gamma or neutron interactions with full charge 
collection [3]. This selection underestimates the actual total population of in- 
complete charge collection events. In order to provide an order-of-magnitude 
estimate of the rate of surface betas, we assume a 50% selection efficiency as 
hinted by Fig. 13 of Ref. [2j. The results are given in Table 2: the count rates 
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are similar to those of alpha particles. Furthermore, the same ratios between 
the counting rates in each detector are observed. We also notice that the in- 
termediate population appears to display energies below ~ 60 keV, which is 
the maximum energy for an electron from a ^^"^Pb beta decay. These obser- 
vations point toward a correlation between the identified alpha population and 
this intermediate population. 

4.4 Contamination scenario 

In conclusion, we have identified three event populations consistent with a single 
coherent contamination scenario: 

• A population of alpha interaction events from the disintegration of ^-'^"Po 
very close to the detectors, with rates of the order of 5 counts/kg/d in the 
fiducial volumes of the detectors. 

• A population of ionization-less events, with energies below 100 keV, which 
we can identify with ^''^Pb recoils originating from the copper surfaces 
facing the detectors. The count rates are of the same order as those 
observed for alpha interactions, and the rate ratios between detectors are 
consistent with alpha data. 

• A population of events with ionization yields lying between those for 
gamma rays and those for nuclear recoils, associated to surface interac- 
tions, which we link to electrons emitted in the decay of ^^''Pb. Here again, 
the rates are of the same order as those observed for alpha interactions, 
and the ratios between detectors are conserved. 

This leads us to strongly suspect the ^^°Pb contamination of the copper parts 
facing the detectors while the hypothesis of a contamination of our detectors 
themselves is not favoured due to the absence of registered ionization-less events 
with a nominal energy of 103 keV. 

It is worth mentioning that the gamma component of the background pro- 
duces also Compton electrons which can escape from the copper surfaces sur- 
rounding the detectors. MiscoUected surface events can be produced when such 
low energy electrons hit a detector. The simulation of the dominant gamma 
background (from the bulk of the copper shield) has shown that, below 200 keV, 
0.6 % (1.1 %) of the events are electrons interacting in the first 20 /im (50 /im) 
under the electrodes; electrons coming from outside of the detectors consti- 
tute 70 % of these populations. The absolute rates are between 2 counts/kg/d 
(20 ^m) and 4 counts/kg/d (50 /im), the actual fraction of miscoUected electrons 
depending on the variation of the charge collection efficiency with depth. 

Finally, a complete model of beta contamination on the detector and copper 
surfaces should also consider ^^C. All surfaces are usually quickly covered by a 
thin (about 1.5-4 nm) layer of organic compounds CHj, (see for instance |17j). 
Given an isotopic ratio ^'*C/^^C= 1.3 x 10^^^, this yields 0.5 to 3 beta elec- 
trons/kg/d for the central part of a detector, with a 156 keV end-point. Com- 
parison with "intermediate event" rates given in Table 2 shows that there might 
be a contribution from ^"^C contamination. 
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5 Neutron background 



Fast (0.5-10 MeV) neutron interactions are a dangerous background as they 
produce nuclear recoils that can mimic a WIMP signal. One coincidence between 
nuclear recoils in two detectors due to a neutron was observed in EDELWEISS-I 
data [2] . We present in this section the simulations of this neutron background. 

5.1 Simulation of neutron transport in EDELWEISS-I 

Before any attempt to evaluate the number of nuclear recoils in the EDEL- 
WEISS - I data due to interactions from a residual neutron flux, it is necessary 
to test the reliability of the Monte Carlo code used for neutron transport through 
a simulation of a calibration run with a ^^^Cf source. Using GEANT3 a nor- 
malization factor of 0.68 had to be applied to the simulation to reproduce the 
experimental rates [18j . An alternative simulation was then performed with 
MCNPX [12]. The detector response was applied to the simulated germanium 
recoils as described in Ref. [S] . The data and the simulation were then analyzed 
using the same selection procedure to extract the recoil energy spectra. 

The ^^^Cf source was positioned, through a 50 mm diameter hole in the 
lead/copper shield, at 70 mm from the external wall of the lead shield and at 
20mm below the detector stack: the neutron count rate decreased from the 
bottom (GGA3) to the top (GSA3) detector. These rates are well reproduced 
by the simulations as can be seen in Fig. [7] The ratio MCNPX/data of total 
number of events is 1.067 ± 0.079 (measured activity of the source [IH|) ± 0.015 
(Stat) ± 0.007 (MC Stat). The measured proportion of events in the 3 detectors 
of 27%, 32% and 40% is well reproduced by the simulation (30%, 32% and 40%). 
However the shapes of the energy spectra differ slightly, a feature also observed 
with GEANT3 simulations. 

5.2 Neutron flux in the Modane underground laboratory 

With a rock overburden equivalent to 4800 m of water, the residual neutron 
background originates mainly from the radioactivity of the rock. Two contri- 
butions arise from the contamination of U/Th: spontaneous fission and (a,n) 
reactions. In Ref. [20], the neutron flux coming out of the LSM rock had been 
deduced from the measurements made with a ^Li-doped liquid scintillator cell, 
using GEANT3 for the neutron transport simulation. This spectrum was in 
turn interpreted in terms of the sum of the contributions from spontaneous fis- 
sion and (a,n) reactions, using a simplified calculation for the latter process. In 
Ref. [21], the original data had been re- interpreted in the light of more reliable 
neutron transportation codes (optimized versions of GEANT3 and MCNPX), 
having for consequence the reduction of the estimated neutron flux in the LSM 
from 4.0 10~^ to 1.6 10~® n/cm^/s. In the present work (flrst introduced in 
Ref. i22i), we conclude these studies by using the SOURCES ^ code, which 
includes a more exact calculation of the (a,n) contribution and using MCNPX 
for neutron transport. 

The neutron flux in the LSM was simulated using the procedure described 
in [21] . The SOURCES code was used to calculate the rate and energy spectrum 
of neutrons produced by spontaneous flssion and (a,n) reactions due to the 
U/Th contamination in the rock. The LSM rock composition and its U/Th 
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contamination was taken from Ref. |20| . The neutron propagation in the rock 
to the LSM cavity takes into account backscattering on the walls. In the present 
calculation, the concrete covering the walls is not present. 

The measurement of Ref. (50] was performed with a ^Li-doped NE320 liquid 
scintillator cell, where neutrons were positively identified from the observation 
of a proton recoil followed by the neutron absorption on ^Li. The detector 
(8.5x8.5x85 cm^^ scintillator) and its shielding (5 cm Cu and 12 cm Pb) were 
simulated with the MCNPX code. The detection efficiency is about 10% for 
incident neutrons above an energy threshold of 2 MeV. The following treatments 
were applied: light yield efficiency based on Birks' law fitted to proton data, 
10% energy resolution, 43% event selection efficiency [20]. 

Fig. [5] shows the comparison of the simulations to the data. The SOURCES 
spectrum agrees better in shape with the data than the spectrum determined in 
Ref. [10]. Note however that all simulations are normalized to data. A normal- 
ization factor of 2.28 is needed in the simulation, which could be explained by 
inhomogeneities in the rock composition of the LSM walls, especially its water 
content. The large influence over the neutron flux of a small amount of hydrogen 
(mainly in water) in the LSM rock is emphasized in Ref. ,24] : the 1% hydrogen 
content reduces the neutron flux above 1 MeV by a factor 2.1 (see also Ref. dS] 
for a discussion of this effect in the Gran Sasso context). 

Fig. [9] shows the results of the simulation for the neutron spectrum in 
the LSM cavity after normalization to the data. We take the normalization 
factor as a systematic uncertainty on the simulation. We obtain a flux of 
1.06 ±0.10(stat.) ±0.59(syst.) 10"^ n/cm^/s above 1 MeV. Assuming a fuU cov- 
erage of the walls by concrete and its U/Th contaminations taken from Ref. [20] 
would give essentially the same result on the flux (less than 1% higher) with a 
corresponding normalization factor of 2.03. 

5.3 Neutrons in EDELWEISS-I 

The rate of nuclear recoils due to the neutron flux coming out of the rock, as 
estimated in Fig.jHl has been calculated using MCNPX. This flux is transported 
through the 30 cm paraffin shielding and the experimental setup shown in Fig. 1. 

We expect about 0.026 ±0.002(MC stat) ±0.018(syst) ncutrons/kg/d from 
the rock radioactivity, 0.002 ±1.2% (MC stat) neutrons/kg/d from the 0.25 ppb 
^^*U contamination of copper shield and less than 0.001 neutrons/kg/d from the 
upper 0.1 ppb limit on ^ssy contamination in lead shield. This translates into 
about 1.6 ±0.1 (MC stat) ±l.l(syst) nuclear recoils expected in EDELWEISS - I 
data (62 kg-day). Fig. fTUl shows the corresponding recoil spectra in the detec- 
tors. The experimental spectrum of Ref. [2], with 34 events between 15 and 
30 keV and 3 betwen 30 and 100 keV, can easily accommodate the presence 
of a few nuclear recoils due to neutron scattering. The only direct experimen- 
tal proof of the presence of a neutron flux is the observation of a coincidence 
between nuclear recoils in two detectors [2]. The present simulation gives fur- 
ther support for this interpretation. First, the rock spectrum of Fig. fTUl shows 
that the recoil energies in this coincidence (14.8 and 14.5 keV) are typical for 
neutron scattering. Secondly, the simulation confirms that the ratio of single 
to coincidence is approximately 10:1, consistent with the ratio observed in neu- 
tron calibrations. With this ratio, it was concluded in Ref. [^ that from one 
coincidence, the prediction range for accompanying single events is from zero to 
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40 events at 90% CL. Although consistent with the hypothesis of a baekground 
originating from both surface events and neutron interactions, this range does 
not help constrain the relative importance of these two contributions. 

6 Conclusions and prospects 

Three main backgrounds have been identified in the EDELWEISS-I data sets. 
The sensitivity is limited by background events which, after the nuclear recoil 
selection, mimic true WIMP induced nuclear recoils with a low value, about 
0.3, of the ionization to recoil signals ratio Q. EDELWEISS-II, the second 
phase of the experiment, is designed for a two orders of magnitude sensitivity 
improvement relying upon the efhcient rejection of these background events. 

Before the rejection of the bulk electron recoils, the overwhelming majority 
of the events is a gamma population entirely dominated by the U/Th contam- 
ination of the Cu shield. These events are rejected with an efhciency greater 
than 99.9%. Nevertheless Compton electrons can escape from the closest copper 
surfaces, reach the detectors and possibly produce miscoUected surface events; 
the gamma component has then to be maintained as low as possible. In the 
EDELWEISS-II set up, extensive material selection, mounting and operation of 
the cryostat under clean room conditions (class 100), secure a better radioactive 
cleanliness. Copper is no longer used for shielding. The inner part of the shield 
is made of very low radioactivity archeological lead. Based on GEANT3 sim- 
ulations and the results of the activity measurements, the overall gamma rate 
for this new setup is predicted to be between 1.0 and 0.1 of the gamma rate of 
EDELWEISS-I, using upper limits or central values, respectively. 

A very small ^^°Pb contamination at the surface of the copper detector cas- 
ings is the very likely source of the observed alphas and near-electrode electron 
events showing a deficit of the charge collection. This probably comes from 
an exposure to radon at some steps of the manufacturing and handling pro- 
cess. In order to minimize the radon exposure, the EDELWEISS-II clean room 
is supplied, during the detector mounting phase, with air of very low radon 
concentration delivered by the LSM radon-trap facility. Nevertheless a more 
decisive approach consists in the identification and rejection of the near-surface 
events. Detectors equipped with Nba;Sii_a; thin films as thermal sensors have 
been operated in the last months of the EDELWEISS-I phase. These films are 
sensitive to the transitory athermal part of the phonon signal, which consti- 
tutes a near-surface tag [26j. These first tests with 200 g modules have shown 
a factor of ten reduction of the surface event rate while retaining a 80% effi- 
ciency. Seven 400g NbSi modules, over a total of 28 detectors, will be operated 
in the first phase of EDELWEISS-II. Other possible solutions are still in a R&D 
phase: identification of surface events using interdigitized electrodes or pulse 
shape analysis of the charge signal . 

The last identified background component comes from neutrons, which, as 
WIMPs, induce nuclear recoils through elastic scattering. The corresponding 
count rate is not yet much constrained by EDELWEISS-I, but its existence 
is established by the observation of one coincidence event. For EDELWEISS- 
II the reduction of this neutron background becomes critical. A 50 cm thick 
polyethylene shield all around the experiment moderates the low energy neu- 
trons. Neutrons created by muon interactions in the shielding are tagged with 
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a 5 cm thick plastic scintillator muon veto of 100 surrounding the whole 
experiment (95% coverage). The compact arrangement of the multi (up to one 
hundred) detector structure allows further rejection through anti-coincidence 
between detectors. 

EDELWEISS-II is now running at the Laboratoire Souterrain de Modane 
and first results are scheduled for 2007. 
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u 


Th 


4Uj^ 


'^'^Co 






(ppb) 


(ppb) 


(niBq/kg) 


(niBq/kg) 


(mBq/kg) 


Copper shield 


0.25 ±0.06 


0.44 ±0.27 


< 15 


< 0.6 


300 ± 150 


Copper holders 


< 0.1 


< 0.1 




< 1.0 




Cubl springs 


< 0.6 


< 1.7 


40 ±30 




260 ± 130 


Ronian lead 


< 0.022 


< 0.032 


< 1.3 




< 200 


Teflon balls 


< 0.2 


< 0.5 


80 ±30 


5±2 


40 ±20 


Wires 






1400 ± 1000 







Table 1: Measured contaminations or activities for several materials used in 
EDELWEISS-I. Statistical errors are at 1 cr and limits at 90% confidence level 
(the detector background at 46.5 keV leads to high errors on ■^^"Pb activities). 



Detector 
Electrode 


GSA3 
Center Guard 


GSAl 
Center Guard 


GGA3 
Center Guard 


Alphas count rate (/kg/d) 


5.0±0.8 24.7±2.2 


5.2±0.8 17.8±1.7 


2.4±0.6 13.3±1.5 


Heavy nuclear recoils 
Count rate (/kg/d) 


5.4 ±0.8 


2.1 ±0.6 


1.5 ±0.5 


"Intermediate events" 
Count rate (/kg/d) 


6.3±1.0 33.0±2.5 


5.0±0.8 33.3±2.4 


3.2±0.7 20.6±1.9 



Table 2: Count rates for alpha events, heavy nuclear recoil events (recoil energy 
greater than 40 keV) and "intermediate events" (between 15 and 200 keV recoil 
energy). These latter are defined as the events between the electron recoil band 
at 3.29cr and the nuclear recoil band at 1.65(7 (see Fig. 6). 



Particle 


Energy 


Cu 


Ge 


Pb 




10 keV 


9 /im 


170 /im 


18 /tm 


Gamma 


100 keV 


6 mm 


8 mm 


400 /im 




1 MeV 


40 mm 


80 mm 


30 mm 




10 keV 


200 nm 


350 nm 




Electron 


100 kcV 


11 /im 


20 /im 






1 MeV 


340 /im 


700 /im 




Alpha 


5.3 MeV 


11 /tm 


19 /im 


15 /im 


Polonium 


100 keV 


40 nm 


68 nm 





Table 3: Typical penetration lengths of various particles in Cu, Ge and Pb. 
Values for gammas correspond to a 10 % transmission probability. Values for 
electrons correspond to the maximum depth at which an electron has deposited 
90 % of its energy (results from CASINO \l2i simulations). Values for alphas 
and Po nuclei are the mean penetration length given by SRIM [T3| simulations 
with a normal incidence. 
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Figure 1: Schematic view of the EDELWEISS-I cryostat within its Cu and Pb 
shields, as implemented in GEANT3 Monte Carlo simulations. Orange areas 
figure copper, grey textured areas stand for lead. The inset shows the three 
germanium detectors (from top to bottom: GSA3, GSAl and GGA3) encased 
into individual copper casings. Teflon balls and Cu springs are visible in the 
upper right corner of the detectors holders. Wires to the cold electronics are 
going up at the left of the detector stack. Not represented on the figure is the 
30 cm external paraffin shield against neutrons. 
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Figure 2: Schematic cut of an EDELWEISS heat and ionization germanium 
detector. The thickness of Al electrodes, amorphous layer and NTD sensor are 
not represented to scale. 
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Figure 3: Gamma background energy spectrum (ionization channel) compared 
to Monte Carlo simulations of various sources of radioactive contamination. 

(a) The two experimental data sets 2003h (full line) and 2003p (dashed line). 

(b) Sum of the 2003p and 2003h data sets (full line), simulations of the U/Th 
(dashed hue) and ^'^Co (dotted hue) contents of copper close to the detectors (de- 
tector holders, cryostat). The quoted radioactive contaminations are measured 
upper hmits. (c) Sum of the 2003p and 2003h data sets (full hne), simulated 
contribution of the radon trapped in the lead-copper shield (upper limit, dotted 
line), (d) Sum of the 2003p and 2003h data sets (full line), simulation of the 
U/Th content of the copper shield (measured concentrations, dashed line). 
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Figure 4: Distribution of the ionization quenching ratio Q with recoil energy Er, 
for the data of the run 2003h. Bottom left: Ej^ projection for 0.2 < Q < 0.4. 
Bottom right: Q projection for 4.5 < En < 6 MeV. 
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Figure 5: Amplitude distributions of ionization-less events after rejection of 
NTD events and noise events. The heavy nuclear recoil events at the surfaces 
belong to this population. 
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recoil energy (keV) 

Figure 6: Representative sample of the "intermediate event" population: run 
2003p, sum of the fiducial volumes of the three detectors, 22.5 kg.d. Full (red) 
circles: selected "intermediate events" ; (black) dots: remaining part of the low 
background data. 
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Figure 7: Energy spectra of the nuclear recoils in EDELWEISS-I calibration 
runs with a ^^^Ci neutron source. The MCNPX simulation is also shown (full 
line, no normalization to data). 
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Figure 8: Neutron flux measurements in 
the LSM cavity from Ref. [55]. The 
data points are the experimental electron 
equivalent energy spectrum. The dot- 
ted line is the original simulated spec- 
trum from rock radioactivity [20j . The 
full line corresponds to the present work 
using SOURCES and MCNPX. Both sim- 
ulations are normalized to the data. 
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Figure 9: Simulated neutron energy spec- 
trum in the LSM after normalization to 
data of Ref. [2D] . The neutron production 
from 238^ g^^^ 232rp]^ traces in the rock 
is simulated with SOURCES and propa- 
gated with MCNPX. 
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Figure 10: Expected energy spectra of the nuclear recoils in EDELWEISS-I low- 
background runs from radioactivity of the rock (full line), ^ssy contamination 
in copper (dashed line) and in lead (dotted line). 
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